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SUMMARY 

A new method is described for estimating low speed wind tunnel blockage, 
including model volume, bubble separation and viscous wake effects. A 

tunnel-centerline, source/sink distribution is derived from measured wall 

pressure signatures using fast algorithms to solve the inverse problem in 

three dimensions. Blockage may then be computed throughout the test volume. 

The method is applied to the results of tests carried out on normal flat 
plates, on spheres, on an idealized automobile, on flat plate wings, and on a 
powered aircraft model. Correlations using scaled models or tests in two 
tunnels were made in all cases. In many cases model reference area exceeded 
10% of the tunnel cross-sectional area. Good correlations were obtained 
regarding model surface pressures, lift drag and pitching moment. 

it is shown that blockage- i nduced velocity variations across the test 
section are relatively unimportant but axial gradients should be considered 
when model! size Is determined. 

Detailed descriptions are given of the computational procedures and a 
program listing is provided. 

A copy of this document is retained in the Lockheed -Georgia Company 
Engineering Report Files. The identifying number is LG79ER0051. 


' klceixn'g page blank not filmed 


. »f I 


XI I 


ACKNOWLEDGEMENTS 


The authors wish to offer their thanks and apprec I a t i on to the follow! 
people for their help in this program: 

Vic Lyman for calculating the influence coefficients for the AO' x80' 
tunnel cross section. 

“Sam" Sampath for studies of sphere aerodynamics using the Douglas- 
Neumann program. 

Sherry Bartley for help in preparing the figures. 

Barbara Reagan and Elizabeth Gardner for typing the test and figure 
mater ia I . 


1.0 INTRODUCTION 


1.1 Reasons for the Present Study 

During the course of routine testing in the Lockheed-Georg i a Company’s 
l6i’ wind tunnel, the shortcomings of existing wind tunnel blockage 

correction methods became evident. The number of corrections applied roughly 
equalled the number of customers requesting them and when applied to the 
same data, a very significant spread in results was apparent. A review of the 
methods themselves revealed that they were often based on doubtful assumptions 
or requ i red test information which was not easy to obtain. A long-term 
program was therefore started with the aim of overcoming these deficiencies. 

A nev/ blockage estimation method v/as developed in 197^~5 (ref. l) based 
upon the analysis of wind tunnel wall pressure signatures. The method permits 
the estimation of both solid blockage (including that due to closed-separation 
bubbles) and wake blockage for three-dimensional models. In its early form 
"desk-top”, chart look-up techniques were employed and a number of somewhat 
arbitrary assumptions were required. Nevertheless, good correlations were 
obtained in tests with large square flat plates mounted normally to the air- 
stream, Over the intervening years, the method has been refined somewhat 
(refs. 2 and 3) and has been computerized in an efficient form, it is no\^ 
used regularly at Lockheed-Georg ia both in its I6ix23i' low-speed wind tunnel 
and in a small, 30" x A3" research tunnel. The method is proving invaluable in 
full-scale car tests, in the large tunnel, and is used routinely by major 
car manuf o*:.turer s . 

Though the early normal flat plate tests underpin applications to cars 
quite well, a much broader correlation base is required before the method can 
be considered for general use. To this end, a number of basic tests were 
conducted in 1975~6) supported both by Lockheed and by the NASA, involving 
more complicated configurations. A major objective of the present study is 
to analyze and correlate these data and to investigate practical problems and 
limitations which are found. An equally important objective is to prepare a 
program for the NA3A-Am-.s 40 ‘ x80* tunnel, similar to the program now In use 
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.t Lockhe.d and implement It. off line, on the CDC 7600 at A„s, This program 
mented heiein. These objectives relate predominantly to defining 
the model environment, in terms of tonne! - induced velocity increments doe to 
blockage effects. The original terms of reference also emphasised the use of 
small-model assumptions and effective node] position. However, this has been 
found .mpractical In some cases, since large models are needed to produce 
tunnel effects which can be measured accurately. 

1.2 Physcial Limitations for Testing 

Tunnel blockage, due to the volume of a mode! and its wake, causes 
velocity changes orimarily in the axial direction. These changes are rela- 
tively uniform at a given transverse plane (see Reterence 3) so transverse 
velocities are small. However, axial gradients of axial velocity can be 
significant in relation to model length and can lead to increasing uncertainty 
•n correcting data for large models. Adverse tunne 1 - 1 nduced pressure gradients 
may cause noticeable boundary layer thickening or even changes in separation 
point on the model. Even though the new, wall pressue signature methods for 
blockage estimation provide a good definition of interference velocities and 
even though these can be translated into "horizontal buoyancy" corrections 
tunnel- induced changes in the viscous effects on the model (particularly 
separation) may prevent valid testing. 


Another consequence of tunnel blockage stems from the pressure distur- 
bances induced at the wall. Generally, there is a favorable pressure gradient 
to a point just aft of the model, followed by an adverse gradient to an 
asymptote determined by wake diameter. Tunnel wall boundary layers thicken in 
this region and the displacement effect adds to blockage. The authors are 
aware of no previous attempt to correct for this effect. However, the wall 
pressure signature method provides this correction automatically. 

Though, in principle, models with very high solid/bubble blockage and 
low wake blockage might separate the flow on tunnel surfaces, this has not been 
observed in the unpowered tests described herein. Even with a 15- inch super- 
critical sphere in a 30-inch high tunnel (13-5% S/C), there was no sign of 
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tunnel floor or roof separation. It appears that lift effects, particularly 
with powered nrode 1 s , are more likely to separate tunnel surfaces than Is 
b 1 ockage . 

Model lift or any net force transverse to the flow, causes tunnel- 
induced changes at the model location primarily in flow angle. Though the 
topic Is outside the scope of the present work, transverse gradients in this 
flow angle may provide an earlier test limit than do blockage effects. A 
large span wing, for example, may develop un represen ta t i ve span loads, in- 
correct induced drag and even premature tip stall long before a limiting 
blockage condition is reached. 

Whether lift or blockage is limiting, it is unlikely that hard, generally 
applicable test boundaries can be drawn, “Hard points'* can be established for 
represen tat i ve con f i gurat I ons , as herein, but in novel or marginal circum- 
stances direct observation should always be used to determine IP the tunnel 
is “in trouble." As the relationship between compu ta t I ona I and experimental 
approaches to fluid problems become more developed, it will be possible to use 
the new, wall pressure signature methods with increasing e f fee t i venes s . The 
me r ge r v^)il) probably start at a basic re sea rch I eve 1, in applicatt on to si mp 1 e r 
shapes and will spread more widely as the art develops, 

1.3 Layout of the Present Report 

The wo^k presented in this report concerns theory, its application to 
practical problems at several levels, and ensuing correlations of pressures 
and forces, A brief overview of the method and the tests Is given in Sections 
2 and 3 followed, in successive sections, by accounts of successive steps in 
reduction and analysis. At each stage, appropriate examples will be drawn 
f rom the body of test data. Because of the breadth and depth of the present 
subject matter, no attempt will be made to report the data comprehensively or 
to provide extensive theoretical development where these are available 
e 1 sewhere. 
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2.0 OUTLINE OF THE "WALL PRESSURE SIGNATURE" METHOD 
AND ITS RELATIONSHIP TO OTHER METHODS 


2.1 Determination of Solid and Wake-B i ockage Parameters 

In a typical experiment (Figure 2.1) wail pressures are mea jred along 
the centerlines of the tunnel walls or roof. If lift is present, the mean of 
corresponding floor and roof pressures can be taken, or sidewall pressures may 
be used. These pressures are used to determine source and sink strengths, 
spans, and locations on the tunnel centerline (Figure 2.2) which define a body 
outline which is equivalent to the test model and its wake. 


A typical wall supervelocity signature (Figure 2.2, lower) includes a 
velocity peak, just aft of the model, which reflects the presence of the ex- 
panding separation bubble which closes further downstream leaving a viscous 

wake. The asymptote downstream is a result of the displacement thickness of 
this wake. 


For the purpose of ana'ysis, the measured profile is resolved into 
symmetric and ant i -symme t r i c parts, as shown in the lower part of Figure 2.2. 
This simplifies the task of determining the equivalent body since the symmetric 
part coi responds to solid/bubble blockage and the anti-symmetric part reflects 
wake blockage. Resolution of these components from the measured signature is 
achieved via an iterative process which will be described in Section k. 

The wake blockage parameters, and X 2 , emerge directly from the wake 
signature analysis. However, to obtain the solid blockage parameters Q and 
c^, a chart look-up technique is employed as illustrated in Figure 2.3. Prepa- 
ration of these charts is discussed in Appendices IV and V. Inputs to the 
charts are the symmetric part peak height (Au/U) and half width at half peak 
height, Ax/B, from the resolved signature, together with source/sink span 
(b^/B) which is estimated from model geometry. 

In essence, the above analysis is an engineering solution of a set of 
equations involving five nonlinear terms, (source and sink geometry) and two 
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linear terms (source and sink strengths). These equations are written out 
fully, for rectangular tunnels, in Appendix IV. 


Because the problem is nonlinear, multiple solutions are possible and it 
is necessary to demonstrate that appropriate roots are obtained by the above 
procedures. Reference 3 describes a Lockheed study using a nonlinear solver 
algorithm which accepts sets of seven wall pressures as inputs and returns 
the seven aerodynamic variables as output. It was necessary to provide the 
program with an initial guess of the roots. By varying the guesses, it was 
possible to converge onto various differing roots. By these means, it was 
established that the source/sink strengths were "hard" variables (i.e., had 
nearly-unique roots) and their streamwise locations were "soft" in nature, 
table at the top of Figure 2.4 shows a number of sets of roots obtained for a 
typical normal-plate experiment. The lower part of the figure shows the 
corresponding centerline distributions of interference velocity. Despite the 
wide variation in the aerodynamic variables shown in the table, the plot shows 
that this was inconsequential with regard to the interference velocity at the 
tunnel centerline. The engineering method produces solutions very close to 
those from the full nonlinear analysis. 


The 


2.2 Interference Velocities and Their Distribution in the Tunnel 

Figure 2.5 shows incremental pressure coefficients, at four tunnel cross 
sections, Induced by the image system for a large sphere at subcritical 
Reynolds number. The corresponding supervelocities, on the tunnel centerline, 
increase from 7-5% at the front of the sphere to 3.8% at its maximum diameter 
and continue r*- ‘ng to 11.3% and 11.2% at planes C and D due to the expanding 
separation - >ble. Interference velocity decreases downstream of D as the 
sepa'^atton bubble closes. 


..lOugh axial variation is very marked, the lower part of Figure 2.5 
shows that percentage variations at any one cross section are not large - 
particularly within the region occupied by the model. It follows from 
conUnuity considerations that b lockage- i nduced lateral velocities - which 
derive from axial velocity changes - are also small. A detailed analysis of 


blockage variations on sphere pressures will be presented in Section 6 of 
this report. 


2.3 Comparison With Other Methods 

Blockage correction methods fall broadly into two classes, requiring 
geometric or geomet r i C“p 1 us-measured tunnel data, respectively. In the ab- 
sence of pressure data from the tunnel, aerodynamic estimates have to be made 
which are often conjectural. Figure 2.6 shows a comparison between five 
blockage estimation methods applied to bluff models, using data taken largely 
from Reference 6. The methods tabulated in the figure will be reviewed only 
briefly here; Reference 6 describes the more subtle distinctions between them. 


It is evident that the two geometric methods (References 7 and 8) give 
lower blockage estimates than two of the three methods employing tunnel 
measurements . Hensel's method, published in I 951 (Reference 9) relies upon 
wall pressure measurement as does the Hacke t t/Wi 1 sden method. However, it has 
been found that use of pressures "opposite to the model location" as recom- 
mended by Hensel may be hazardous, since strong pressure gradients are usually 
present at this location. The method has recently been revised, for use in 
car testing, so as to use peak suction measu remen t s . Another feature of the 
Hensel me .hod is that it uses doublets to model the flow. This renders the 
method less flexible than the Racket t/Wi I sden method, which permits variable 
source-sink spacing. This can be important for long models or models with 
thick wakes. Corrections via the Hensel method (Figure 2.6) are generally 
lower than by Hackett/Wi 1 sden —noticeably so in the case of the sphere. 


Haskell’s method, which appeared in I 963 (Reference 10), determines 
dynamic pressure correction from an empirical relationship involving measured 
drag or, in a variant of the method, measured drag in combination with base 
pressure. Haskell shows the drag coefficient is proportional to base pressure, 
for square plates, and notes that non-square plates exhibit varying base 
pressure across their span. In the equation* Aq/q = nC^ S/C, n is assigned 


xvoid xoifusioyi bloakagc xelocitu vatio, Maskell^s e 

r^eplaaed by n. 


tabulated values which vary Prom 2.77 to 2.13 for plate, aspect ratios in the 
range 1 to 10. Haskell's method is well founded and gives good results for 
normal flat plates. However, in other applications determination of base 
pressure for direct use in the method may be hazardous not only because of 
the aforementioned aspect ratio effects but also due to vortex separations 
which may exist to each side of the separation bubble. If the tabulated 
values of n are used, there may be difficulty both in separating bubble drag 
from other components and in estimating the aspect ratio of the separation 
bubble. Though Haskell asserts correctly that an n value of 2.5 ± 0.25 is 
reasonably accurate in situations where q-correct ions are of order 10°4, the 

spread m estimated Aq/q becomes unacceptable for large models (see Figure 
2.6, bottom ] i ne) . 


It has been noted already that the Hackett/Wi 1 sden 
garded as an extension or generalization of the Hensel me 
shown that the present method, implemented via the "e 
^Reference 3), exact’y parallels Haskell's approach while 
pressure measurements and estimation of g. The relation 
two methods will be examined further as part of the drag 


method may be re- 
thod . I t may also be 
" approach 
avoiding both base 
ship between the 
analysis in Section 7- 


3.0 OUTLINE OF REFERENCED TESTS 


3.1 General Comments 

Figure 3* 1 1 ists tests from which data was extracted for the present 
study. This list Is not a I I - i ncl c s i ve since certain of the tests listed are 
second attempts following earlier d i f f i cu I t ies . For example, all data quoted 
herein for the 30“ x 43" tunnel were taken after the stretched working section 
was installed; the earlier test section length was inadequate. The 1974 flat 
plate tests In the l6i‘ x 23i’ tunnel employed an aft“Sting mount (Figure 
3.2), but were preceded by tests using vertical supports. Some of the latter 
data are quoted in Reference 1. 

For the simpler models, a number of geometric copies were made and 
tested in the same tunnel. The smallest models were sized to give essentially 
free air data and the largest had an area of 10^ or more of the tunnel cross 
section. Tunnel speeds were adjusted to yield the same Reynolds number range 
at each model size. Definitive data such as those quoted here have not been 
found in the literature for large models. 

Powered models are generally too expensive for the multiple-model 
approach to be used. Scaling problems may also be encountered. Accordingly, 
large tunnel versus small tunnel comparisons were used, recognizing that 
differing turbulence cha racter I s t i cs , for example, can cause anomalies. 


3.2 Tests on Unpowered Models 

Mormal flat plates were tested in the l6i' x 23i' tunnel for S/C values 
ranging from 2.38% to 9-53%. All had sharp edges, chamfered from the rear and 
were s t i ng-moun ted centrally from the rear (see Figure 3.2). The horizontal 
sting was mounted to the standard, virtual center tunnel balance via a vertical 
strut which was shielded from the mainstream. Base pressures were measured 
both on the plate aft surface and within the bubble, using a boom (see Figure 
3.2) introduced only during pressure measurements. 
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Ofher normal plate tests are included as a special case of the flat 
plate wing tests (Figure 3.3). These plates had a half-wing aspect ratio of 
1.5 and ranged from 1.56* to lfc.7*«, of the 30" x 1.3" tunnel cross-section. 

An under-fioor balance was again used, great care being taken in scaling the 
base plates and floor clearances. No attempt was made to remove the floor 
boundary layer and, in consequence, data for the smallest model were inferior 
to the remainder. Increased speeds, used to equal I ae Reynolds number across 
the model spectrum, evidently failed to thin the floor boundary layer 
adequately for the smallest model, W1. 


For compatibility with the flat plate tests, a 50% mounting location was 
used and measured pitching moments were re-referenced to quarter-chord during 
data reduct. on. The lift performance of the plates, being sharp edged, was of 
course modest and the addition of leading edge roundings had little effect 
Three spheres, of 8.5, 12.0, and 15-0-inches diameter were fabricated. Figure 
3. shows the 15-inch sphere mounted in the 30- x wind tunnel. The spheres 
were constructed from laminated mahogany and hollowed out for lightness and 
access to instrumentation. Each sphere was fitted with Ik surface pressure 
or.f.ces and was mounted on a 3A-inch diameter strain-gauge balance. The 
sphere was chosen as the solid object model (vis-a-vis flat plate) for 
continued blockage method verification tests because its free-air properties 
are well known and are markedly different at sub- and supercritical Reynolds 
numoers. Furthermore, its potential flow aerodynamics are simple analytically. 

e two larger spheres were tested in the Lockheed 30- x i»3" tunnel under 
blocked conditions. All three spheres were tested at essentially free-air 
conditions in the Lockheed l6i' x 23V tunnel. In some of the latter tests 
static pressures were also measured at a location corresponding to the roof’ 
centerline in the smaller tunnel. These data were used to test pressure 
correction procedures (Section 6). 


3.3 Tests on a Powered Model 


Powered model data quoted in 
winged knee-blown flap model which 
Dimensional details are given in F 


the present report involved a straight- 
had a flap upper surface angle of 76 degrees 
igure 3 . 5 . The wing has a span of 76.2 cm 






(30 ) with tips fitted end e tip cho rd of 1Z.7S cm (5'03'*)» ^ f c ren cc 

area of 973-8 sq. cm (l.048 sq.ft.). Tips-off reference dimensions are used 
throughout this report, as previously (Refs, k and 5). These are: span, 50.8 

cm (20"); chord, 10.16 cm (A"); and area 516.1 sq. cm (0.556 sq. ft.). Full 
dimensional details and a detailed model description may be found in Reference 
2. A general view of the model is given in Figure 3-6. The wing was mounted 
on a fuselage, which housed a three-component strain-gauge balance and air 
bridge. The sting was attached to an incidence quadrant mounted beneath the 
tunnel floor. Similar arrangements were used in the 7' x 10' NASA/AAMRDL 
tunnel. In the smal le>- tunnel, floor boundary layer control was available, 
for use at higher blowing rates, employing either a moving ground or a 
tangential ly-blown ground. Descriptions of there devices and their applica- 
tion may be found in References 2 and S- The same references give a detailed 
amount of calibration and test procedures used for the tests quoted herein. 






4.0 MEASUREMENT AND ANALYSIS OF WALL PRESSURE SIGNATURES 


4.1 Measurement of Wall Pressure Signatures 


Surface pressure orifices are installed in the tunnel 
tunnel roof and possibly along the tunnel floor along their 
lines. Longitudinal spacing should be between 5? and 10? of 
with larger spaclngs towards the downstream end of the test 
recommended model position, where an option exists, is about 
section length from the forward end. A relatively long test 
desirable, preferably about 1.5 times tunnel width. 


s i dewa I I s , In the 
respective center- 
tunnel width, 
section. The 
one-th I rd of test 
section is 


Well-made orifices are required, drilled in metal and free from burrs. 
Individual discs may be inserted in large tunnels. Continuous pressure tapped 
strips may be more convenient in small tunnels. An example of the latter may 
be seen in the tunnel roof in Figure 3.4. Discs are set into the far wall. 
Orifices are also drilled directly into the plexiglass windows in this tunnel. 

^ When run time is not at a premium, a single scanivalve is usually 
sufficient for pressure measurements. However, in a production tunnel or where 
long tube runs would cause undue pressure lags, individual transducers may be 
used. The latter approach is used in the Lockheed l6i' x 23i< tunnel. 


If flow unsteadiness is present, long-time averages should be taken over 
a period of several fluctuation cycles. In setting up and calibrating the 
instrumentation it should be noted that excursions may be as little as .02, 
for low-drag models, or as much as 1.5 for large, bluff objects (see Figure 

4.1). With care and good instrumentation, adequate accuracy can be achieved 
without changing t ransducers . 


A datum 
tunnel empty, 
preferable to 


set of tunnel wall pressure measurements should be 
prior to model installation. Experience has shown 
gather these data at the approximate test value of 


made wi th the 
that it is 

q- 


was on 


figure 4., illustretes a „u«,er of features characteristic of the 
ressur, s,g„atures for relative, y si„p,e „,els. The test concerned 
riat plate wing (Fiaure t 7 ) =roh 

various angles of attack. Since Hni<= 

IS a floor -mounted half-model blort'd^no • r ■ ^ 

tunnel roof pressures. ’ '"fot^ticn .ay he derived solely fro. 

Pressurd^s? T '» 

P sure d.sturhances are guite s.al, and reflect predo.inant ly viscous 
half aT^he''! intensifies and noves 

is uni ; tz ' t : 

of attach a e , h r ' 50-degrees angle 

attack are substantial. For the 80 - d^nH nn-d 

is interrupted at x/B = 0 i; cases, pressure recovery 

wool , fr S ^iiMestino possible wall separation. However 

wool tufts showed no evidence of this. 

figune 4.2 shows the effect of i I ft on .easured sidewall pressure 
ignatures. The pressure coefficients ate uncotrected for e.pty tunned 
irregularities, which show most clearly 'n the « n 

height orifirec . the a = 0 signatures. As the mid- 

Ssiailer th ’• ""‘i"" increments due to lift were scewhat 

sialler than ,n the lower corners of the tunnel, nonetheless, the signatures 

own couid be used to estimate „ft interference using n^thods sl.ilaTo 
those described herein for blockage estimation. 

Figure 4.3 provides a further example of the way in which a wall 

pTir^ --- ^ -e'. this time a 

At suocritical Reynolds numbers there is a substantial separation 

nulbe^ th! fLr'irT " itical Reynolds 

This •’ ^ ^ decrease in wall suction is evident 

^ signature is close to that calculated for potential flow. The automati'c 
esponse of the method to such flow changes is a major virtue. 


^•2 Signature Analysis 

The overall procedure for determining blockage parameters is described 
■P -oad terms, in Section 2 and is illustrated in figure 2.,. The p„;l 
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of this sub-section is to describe the methods used in more detail. Figures 
A1 and A2 give an outline of the major elements of the signature analysis 
program. Appendix I gives a description from a user's point of view. 

Generation of Basic !xu/U ^ vs. x/B Signature 

Pressure coefficients measured at tunnel surfaces with the tunnel empty 
are first subtracted from corresponding .mode 1 -presen t values measured at 
approximately the same 'q'. The resulting increment, C . is then converted 
to a velocity increment using the equation: 


Au 


(I - AC )^ - I . 
P 


(4.1) 


Figure 4.4(a) shows a typical measured signature at this stage. 

Resoouvton of Szgnatuve into Symmetric and Antisymmetric Parts 

Resolution into symmetric and antisymmetric parts is an iterative pro- 
cess which starts with a trial estimate of the antisymmetric part. This is 
subtracted from the symmetric part and a test is made to determine whether the 
inflection point of the antisymmetric component coincides with the peak of the 
symmetric component. If the separation is too great, the wake source is moved 
to the peak location and the process is repeated. Up to five iterations are 
permitted. The individual steps in this process are described below. 


Calculation of tke Antisymmetric Signature 

For an ideal signature measured in a long test section, Au/U^ asymptotes 
to zero upstream and to a constant value, determined by wake diameter, far 
downstream. The downstream asymptotic value (Au/Uj,^,^ is noted and an esti- 
mate is made for wake source position x . The model position is used initial- 
ly. These are substituted, together with the assumed wake source span, b /B, 
into the equations ’ 


|4 


where 


/Auv 

^11 ) - A5 (i 




CO A 


5(1 + tanh (Ag (~) + A7)) 


0.5 (Au/U ) 


MIN 


^^. 2 ) 


^6 = 3-30 + 0 . 2 h 7 h EXP (7.8131 b^/B) | 

A? = -Ag (x 2/B) j 

‘■'presents the downstream asymptote of the signature. 

Equations (U.Z) are a curve-fitted approximation to t ,e true li,. -o 
equat, on which is too unwieidy for convenient man ipu i a t ion. The constaltr" 
■hvoived the equation are specific to a particular tunne, sh 
.aloes quoted ahove are for a rectangular tunne, of 0.633 aspect ratio. 

On the majority of occasions, the test section win be too short in 
relation to the model and its wake for rhfo 

ar horh a . asymptote to be well established 

oth ends. Designated points may therefore be selected at the 

and downstream ends of the signature. On the first pas - rh H 

( . -.these are assumed to be asymptotes. "Hiss-d i stancesu are next 

.ecmined, between the designated points and the fitted curve and 

point: "'?hi'; '-P ^«-9nated 

correction is made just once and is not iterated. 

cown " ^> 9 My-powered modeis. If tunnel flow break- 

detects'l'is 

detects this situation by testing the slope of rh« i • ■ • ■ 

downstream point to the preceding point. If .be s I oj i 

P-ticu,ar positive value, a 'walking, procedure is started with , ■ 

until a near-zero slope is found. This establishes rh decreasing 

(A 2) a^ p.n ro 1- estoblishes the asymptote of Equation 

peaks Waha btck? befwaen fha pwo 

biorp, . n " °.e'-eet""ated in such cases and soi id 

blockage is underestimated. 
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When the above process is complete, the fitted antisymmetric signature 
IS subtracted from the measured signature leaving a bell-shaped symmetric 
curve-defined by experimental points [see Figure 4.^i(b)]. 

jf the oyrmeiHe Signature 

A least-squares curve fit procedure is used to fit a parabola to the 
uppermost points of the symmetric signature. Various checks are employed to 
ensure that the points selected both represent the peak properly and extend 
downward sufficiently far chat Ax, the peak half-width, can be estimated with 
confidence. Figure shows a parabola fitted to the upper part of a 

symmetric s i gnature. 

The most significant initial result from the parabolic fit process is 
the peak position x^. If this is sufficiently close to the assumed wake source 
location, X 2 . the source and sink parameters are calculated. Otherwise, the 
wake source location is suitably adjusted and a further iteration is perfoi'med. 

Once the iteration has converged, adjusted values of Ax/B and (Au/U ) 

OO 

are determined for use with Charts I and II. These adjustments correct 
empirically for the fact that the parabola is only a rough fit to the true 
aerodynamic source-sink function. The parabola coefficients are defined by: 

Au/U^ = a (x/B)'^ + b(x/B) + c 

From this we obtain: 


Peak Pos i 1 1 on : 

Xp/B = “ b/2a 



Peak Value: 

(Au/U ) = 1 

max 

.01 (c + 0.5b (x /B)l 
P ' 


Sem i -W j d th : 

Ax/B = 1.045 [ 

-(Au/U ) /2a]^ 

max 

('i.S) 


The factors in {h.k) and (h.S) are the empirical adjustments (see above). 
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Detevminat ion of Flow Model 


Parame tors 


After the stgnature resolution has converged Ax/B and b /B are input Co 
arts I and ’I, via a double interpolation routine, to obtain c /B the 
eource-sink spacing, and (au/ (Q^/B=) which is used to obta in'source 

strength. No„-d i,«ns iona I solid/bubbie source-sink strength is then obtained 

from 


(i., (Au/U ) 

- ^ ^ max 

^^^^2 XIIiZrQ /B2)l 

s • max 


(4.6) 


Source's i'lk positions are given by 


= JP Is 

B B 2 B 


^4 ^ 


B B 2 T 

a/ 


(4.7) 


b /B is assumed to be the same as b^B. which has already been estimated, 
the source sink system is now fully defined. 


SO 


Non-dimensional viscous wake blockage 


source strength is given by 


<1, 


U B2 ~ ^ ^3-^ 

rr> ^ 


where C is tunnel area. 


The wake source position, x^/B and span b / 
signature analysis. 


w 


'B were obtained during 


We now have three source strengths (hQ^/U^B^) , (Q„/U.B=) and (-Q ,'u B=), 

all of s pan ‘ ' ‘ ^ ^ 


tuated at locations (xj/B) . (x./B) and (x^/B) . 
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Detemiyiat ion Contevl 


yie Intevfevonae Signature 


The total centerline signature is built up Prom three components, one 
for each source or sink, as illustrated in Figure ^.5. To determine each 
component at position x/B, the relative source location (x/B - x /3) is 
determined and input to Chart M I to give Au/(Q/B^). Multiplication by the 
appropriate normalized source strength (Q/U^B^) then yields the desired 
interference velocity (Au/U^) on the tunnel centerline. The three components 
are added and the calculation is repeated for each specified x-station. 


Check oyi Curve Fit to 'Wall Measurements 

A recent addition to the program is a routine which determines the entire 
wall signature implied by the flow model parameters determined above. It is 
not essential to determine this curve, but it can be helpful if unfamiliar or 
unexpected results occur. 

The new routine again makes use of the tanh function fit using the 
equa t i on 

Au/U^ = Aitanh (A, | + A3) + Aitanh (-A2 ^ + A4) 

+ Asfl + tanh (A + A )) (i».8) 


where 


Ai = 2 (b2/c) 

A 2 = Ag iStie Equation 4.2) 
A3 = -A2 (xj/B) 
and A4 = +Ao (xu/b) 

A5 , Ag , and A-; are given by equation (k. 2 ). 
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5.0 EFFECTIVE MODEL SHAPE 


5. 1 I nt roduct i on 

For each set of aerodynamic parameters derived from wall pressures there 
IS an implied body shape. It is not necessary to determine outlines of these 
bodies routinely, but useful insight can be gained from them in some cases. 
Calculated shapes can be compared with model dimensions and observed 
separation bubbles. In-tunnel versus free air comparisons can be used to 
make a rough assessment of tunne 1 - i nduced changes in separation bubble shape. 

The equations for the full, two- 1 i ne-source plus line sink bodies are 
complex and most of the present studies have been restricted to the axisym- 
metric shapes generated by point sources and sinks placed in a uniform 
ma instream. 


Continuity considerations lead to the following equation for the body 
radius, r: 


1 3 

f(r) = 1 r2 - y 

Q 

- 

(x - x^) - 

n 

1 + - 

7 L 

n=l 

00 

_ 1 

[(x - x^)2 * r2)i 


where and may be interpreted in terms of ±Q^ and and their positions. 
Thus. Qi and Q 3 become ±Q^ and replaces Q 2 . The x subscripts follow 
accordingly. Equation (5.1) was solved, for r, on an HP29C hand 
calculator using Newton's method in the form 

ff(r. + .) 

r- = _ X I I I 


-1 


where 5. = 10”^ r 

I [ 


The downstream asymptote of the body is given by 


r 

00 


/ ttU 
/ 00 


(5.2) 
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Equation ( 5 .I) is not suitable for determining the location o^ the front 
stagnation point, x . This is given by: 


f(x) = 


3 

' - I 

n= ] 


n 

4ttU 


(x - X ) 

s n 

a¥s (x 

s n ' 


(: 


= 0 


(5.3) 


The above equation is a sextic, having four real and tw. imaginary roots in 
most cases of interest. The lowest x .root corresponds to the front 
stagnation point. The same solver was used as for Equation ( 5 .I). 

5.2 Calculat6d Body Shapes 

For zero solid blockage, equations (5.1 through (5.3) yield the outline 
c: an axisymmetric source body, strength in a mainstream U . Figure 5 1(a) 
shows a number of wake bodies so derived for various wake strengths. Effec- 
t.ve bodies for purely solid blockage (O^ = 0 ) are shown in Figure 5.1(b) 

These two sets of bodies represent extreme cases between which most 
experimental cases will be found. 


F.gure 5-2 shows sample cases which lie between the above extremes. 
Equal values of and produce shapes similar to blunt-ended rods [Figure 

5.2(a)] which develop bulbous ends as relative wake strength diminishes 
[Figures 5.2(b) and (c) ] . 


tach profile in Figures 5 . 


and 5.2 requires nu.nerous solutions to 


Equation (5.1) and this precludes quick, routine estimates, figures 5.3 and 
5.E have ttsrefore been prepared which give the major dimensions of the source- 
source-sink bodies. A reasonable approximation to the body outline may be 
obtained by drawing an ellipse with semi-axes R; and Ax and fairing it into a 
cylinder of radius Rj , obtained from equation ( 5 . 2 ). 


5.3 Comparisons for Tested Models 


Figure 5-5 shows effective shapes 
for the 15 -inch sphere in the 30" x i)3" 


derived from wall p-essure signatures 
wind tunnel. Considering that only 
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three point singularities are used, the effective shapes blend with the model 
outline very well. The change in wake character through transition is well 
depicted and the thinning of the supercritical wake is very marked. It is 
interesting to note that the subcritical wake occupies 2k . 3 % of the tunnel 
cross-section, at maximum diameter, compared with 13 . 7 % for the model itself. 

Generally similar results are obtained for a circular disc mounted 
normally to the flow (Figure 5.6). It is noteworthy that the calculated body 
diameter at the model station closely equals the disc diameter. 

Figure 5.7 shows effective model shapes derived for a less symmetrical 
model, an idealized automobile. It is apparent that when the model is yawed, 
tne present method responds appropriately by shortening the effective body and 
increasing its radius. The ratio of maximum source-sink body area to project- 
ed model frontal area is 0.733 at zero yaw and 0 737 at 30-degrees. If the 
wheel projections are subtracted from the zero-yaw frontal area, ro give an 
aoproximate maximum cross-section, the ratio becomes O.Slk. The supercritical 
sphere had a corresponding area ratio of 0.85k. 

In the subcritical sphere case [Figure 5.5(a)], the calculated wake 
occupies a substantial proportion of the tunnel area and it appears probable 
that interference will distort it. To estimate the magnitude of the change, 
in Equation ( 5 .I) was replaced by a value which included the tunne 1 - i nduced 
axial velocity at the x-station concerned. Figure 5.8 shows that the bubble 
outline is changed surprisingly little. As the model is very large in com- 
parison with those for most routine tests, this suggests that blockage- induced 
distortion of separation bubbles is not likely to be a limiting factor. 

Effects on pressure distribution and implications regarding separation 
location will be discussed in Sections 6 and 7. 
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6.0 CORRECTIONS TO MEASURED PRESSURES 


6.1 Simultaneous, Whole-Model Corrections 

In the preceding sections, methods have been described for defining 
tunnel -induced changes in model environment in terms cf changes of velocities 
u, V and w at points distributed over the surface of the model. As we have 
5een, the predominant effect of blockage is to increase u. A simple increase 
presents no correction problem. Pressure coefficients can be renormalized 
simply and a correction to Reynolds number will ensure that viscous phenomena 
are properly defined. However, there are usually significant changes in 
interference velocity along the model length. In situations involving 
sensitive viscous effects - incipient separation, rapidly moving separation or 
transition - these gradients may render the tunnel data uncorrectable. An 
implicit assumption, in correcting in-tunnel pressure measurements, is there- 
fore that there are no sensitive viscous phenomena and, moreover, that tunnel 
gradient-induced changes in boundary layer thickness are small. 

Within the above framework, there are several levels at which corrections 
may be made. These range from a mathematically represented model surface with 
full and simultaneous application of the boundary conditions to simple super- 
position methods. The former, whole-model corrections will be reviewed in 
this section. Simpler, more practical methods are discussed in Section 6.2. 

To illustrate the application of blockage corrections simultaneously to 
a complete model, a Dougl as-Neuman , source-panel model was set up of the 
15 -inch sphere. The total spherical surface was represented using 6i.8 10 - 
degree by lO-degree panels. For economy, only one quarter sphere (seen from 
the front) was modeled explicitly and imaging techniques were used to model 
the effects of the remaining panels. 


The objective of the ensuing calculation was to calculate in-tunnel 
pressures on the sphere surface and compare them with aieasurement . The source 
and sink image system determined From an experiment at supercritical Reynolds 
number was used to calculate the three components of interference velocity at 
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.ac.pa„„ center T.ase velocities „ec. i„p„t to t.e soetce pane, ptopsa™ 

the rellui'! ' " Figoce 6., shows 

pressures at the sphere equator together with the correspoodip, 

zzzzrT""^""' 

fur he r e " "" '" '*’" -P«F--tal data 

r. et aPt, hut the theor, reproduces the geherai trend weii as far hach as 

the separat.oo point, at about iho-degrees. A sirtiiar correiatio,, was obtain- 

FigureVr""' - - 


An extension of the above study couid be fade in which the effective 

.SCOUS wahe surface [figure 5.5(b), is paneied and blended into the sphere 

Improved correlations in the 120- m i h • spnere. 

tne 120 to 130-degree region would probably result. 

owever, such a detailed approach is not usually justified. 

The above study demonstrates that s imui taneous , whoie model corrections 
red.ct tunnel effects weii for .„deis with predominantly attached flow 

IsL::!: ^ 


6.2 Approximate Method 

Figure 6.2 examines the components of interference velocity and their 
impact upon surface pressures. For this study, a potential flow sphere was 
use . represented by the previous paneled sphere inside the tunnel, and a 
s en, o oublet images to represent the effects of the tunnel boundaries 
IS evident that axial interference components maximize at about 5^ et the 

locaf velocities are greatest - about 0.6“^ - at the 45-degree 

canon. Th, .ntarfarence volocitias ara not tangential to the surface 

r::::r::cr- " ~——f course, „„;„a, 


sets 0 ^ 0 ":: "'"'h “--ponding surface pressure increments. Three 

sets or data are shown, for seoaraf-#^ , . 

w Pirate calculations involving all three inrr^^^- 

-fs (filled circles), transverse components only (triangles) and the axial 
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componen t , 


Au only (crosses). As has already been noted, the 
has a dominant effect. In fact, in a corresponding plot for the 
position, the effect of the transverse components cannot be seen 


Au component 
equator ial 


It is obviously very desirable to correct measured pressures on a point- 
by-point basis, rather than globally, as described previously. Bearing this 
in mind the small spread in points in Figure 6.2(b) and the very large mode] 
concerned, it is also apparent that it will generally be sufficient to concen- 
trate on the An component. Further simplification is afforded if variation 
in Au across the test section is assumed small (Figure 2.5). Only the cen- 
terline increments found in the standard blockage computation are then required. 

The simplest point-by-point Cp correction would be to superpose the 
tangential component of Au on /T^ and recalculate the pressure coeffic- 
ient. This does not work because there is an implicit assumption that the 
mainstream flow component at the sphere surface is affected but the doublet 
strength is not. In the case of the potential sphere, it is clear that, if 
+Au is present at the front of the sphere the stagnation point will move back 
unless the sphere doublet strength increases. It is evident that sphere 
doublet strength must be scaled uo by a factor (U + Au)/U to restore the 
stagnation point to its correct location. Since this factor is also applied 
to mainstream velocity, it is obviously directly applicable to measured sur- 
face velocities, at least for cases with Au constant over the model length. 

Thus, i f q^ is the 'ocal surface velocity in free air, we may write, for 
uniform Au 


in free air and 




CO 




oo 


(1 



oo 


in the tunnel. Equation (6.2) may be written as 


(6.1) 


( 6 . 2 ) 
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Au . 2 
U ^ 


and on substituting into (6.1) this gives 


C 
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c 



(1 




Au v 2 

U ^ 


oo 


+ i . 


( 6 . 3 ) 


I f Au IS not constant along the axis, we may postulate a "local mainstream' 
effect in which Au/U in ( 6 . 3 ) is a function of x. Thus 


Cp (x) 


C (x) - 1 
u 


(1 


~ Au(x) '~2 

I I ' 


+ 1 


( 6 . 4 ) 


To test the method, in-tunnel C^'s calculated with Au, Av and Aw present 
[Figure 6 . 2 (b)) were corrected using Equation ( 6 . 4 ) and compared with a 
corresponding free-air calculation. The results are tabulated below. 


9 ° 

5 

15 

25 

35 

45 

55 

65 

75 

85 

AU/U 

00 

.0385 

.0390 

.0399 

- 0 A ;3 

.0432 

.0454 

.0476 

.0495 

.0506 

*^Pu 

.9775 

.8417 

.5755 

.2208 - 

.1956 - 

.6212 - 

1.0045 - 

1.2955 -I 

.4530 

^Pc 

* .9791 

.8532 

.6106 

.2803 - 

. 1005 - 

.4862 - 

.8303 - 

I.O877 -1 

.2276 

Cp 

Ppree 

.9786 

.8502 

.6034 

.2694 - 

.1133 -. 

.4968 - 

.8350 - 

1.0859 -i 

.2191 

ERROR 

+.0005 

+.0030 

*Using 

+.0126 +, 

equation 

.0109 +. 

4 ) . 

.0128 +, 

,0106 +. 

0047 +. 

,0028 + 

.0085 


The errors in applying equation ( 6 . 4 ) to sphere crest pressures are .013. in 

Cp, at inost. Around the equator, a similar calculation shows errors less than 
.005 at al I posi tions. 
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Bearing in mind the large size of the model used for checkout, it may be 
cone 1 uded that equation (6.4) corrects measured pressures sufficiently 
accurately for most practical purposes. 


6.3 Pressure Correlations on the Model Surface 


Sphere Tests 


Figure 6.3 shows tunnel centerline interference velocities calculated 
from ceiling pres ure d I s tr i but ions for Reynolds numbers through the criticol 
range. At subc ‘ cal Reynolds numbers the tunnel centerline blockage 
velocity increases continuously from about k% to more than 10^ over the chord 
of the sphere. The blockage maximizes just aft of the sphere, at the loca- 
tion of the crest of the separation bubble. The strong variation of blockage 
makes the choice of cha racter i 5 t i c velocity, for Reynolds number determina- 
tion, conjectural. For the present work, the corrected velocity correspond i ng 
to tne suction peak location has been used. A further difficulty arises 
because the subcritical sphere evidently has a significant favorable pressure 
gradient imposed upon it by blockage effects. This raises the possibility 
that the critical Reynolds number will change. 

At supercri tical Reynolds numbers blockage velocities are about half as 
large as in subcritical cases and vary much less over the chord of the sphere. 

Figure 6.4 shows pressures measured on the surface of the 15“ inch sphere 
in the 30" x 43" wind tunnel, designated Cp^ and in the l6i‘ x 23i' wind 
tunnel, designated Cp . Blockage increases both the peak and the wake suc“ 

^00 

tions very significantly in the small tunnel. Correction via equation (6.4) 
leads to good correlation except that large tunnel peak suction is slightly 
higher. However, this is consistent with the slightly higher Reynolds number 
i nvol ved . 

The results of a correspond i ng correction at a supercr i t i ca i Reynolds 
number Is shown in Figure 6.5- In this case the Reynolds numbers were better 
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matched but wake suctions are again slightly higher in the smaller tunnel. 
There are also some detailed differences just aft of the sphere's crest. 

Figures 6.6 and 6.7 show peak suction and base pressure respectively as 
functions of Reynolds number through the critical range. In both cases, the 
corrections to small -tunnel Reynolds number use velocity increments, Au/U^, 
calculated at the peak Cp locations, which are quite close to the sphere 
center in most cases. Pressure coefficients are calculated via equation (6.4) 
using the same value of Au/U^. Application of these corrections increases the 
value of Reynolds number and decreases suction coefficients (see Figure 6.6). 

In both the fully subcritical and the fully supercr i t i ca 1 Reynolds num- 
ber ranges both tunnels indicate similar trends but transition to the super- 
critical character i St i c is earlier in the large tunnel. Hov/ever, examination 
of the base pressure cha racte r i s t i cs , in Figure 6.7, shows no significant 
difference between critical Reynolds numbers in the two tunnels. With the 
exception of the curve for the large tunnel, Figures 6.6 and 6.7 show 

critical Reynolds numbers which agree both between tunnels and with accepted 
turbulence sphere results. 

The above comparisons indicate that for a sphere occupying 13.7^ of the 
tunnel cross-section, application of pressure corrections via equation (6.4) 
produces results comparable with and no less credible than results from 
similar tests in a large tunnel. In correcting both pressure coefficients 
and Reynolds number, interference velocities calculated at the peak suction 
location were employed. Model center values would be a reasonable 
approximat ion to this. 


Test on Idealized kutomohile Shape 

In 1978 , experiments were carried out, under Joint Lockheed/Ford Motor 
Company auspices, on models of an idealized car. Tests were run in the 
Lockheed l6i‘ x 23i' tunnel on a "full size" (frontal area = 5.65?^ of tunnel 
area), on a 0.475, and on a 0.375 scale model (see Ref. 6). Figure 6.8 shows 
the profile in side view. The metal models were very accurately made and were 
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extensively pressure plotted. Design was simplified so that many surfaces 
were flat, being joined by simple radii or sharp edges as appropriate. 

Figure 6,8 shows measured pressures at the centerplane for all three 
models. Floor blowing boundary layer control was used to give a uniform 
velocity profile near the ground and tunnel speed was adjusted to give the 
same Reynolds number in all three cases. Pressure signatures were measured 
along the tunnel ceiling and processed as described previously. 

Figure 6.9 shows the previous data after correction via the present 
method. In this case the model blockage was small enough that a constant 
value of Au/U could be used at all locations. The peak Au was used. Slight 
overcorrect i on is evident in some regions, but the overall correlation 
between models is excellent. 


6.4 Pressure Correlations at Off-Body Points 

In the immediate vicinity of the model it is appropriate to correct 
pressures in the manner just described. However, in the far field, there is 
another option by which off-body pressures may be deduced directly from the 
source-sink system derived from wall pressure signatures To demonstrate 
this, static pressures were measured in the large wind tunnel, in the presence 
of the 15“inch sphere, at locations corresponding to the small tunnel wall 
and ceiling orifice rows. Attempts were then made to predict the large tunnel 
("free air") pressures from the small tunnel data. 

To apply the direct method, source and sink strengths and locations are 
derived from wall pressure signatures, but no attempt is made to correct the 
signatures themselves. Instead, velocities Induced in free air by the 
source-source-sink system are calculated at the desired, off-body, far field 
locations. These are converted to pressure coefficients and give the results 
indicated by crosses in Figures 6.10 and 6.11. At the wall position, correla- 
tion between large tunnel measurements (circles) and the predicted pressures 
is very good. Suction is underpredicted slightly at the ceiling. 
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An alternative method of determining free-air far field pressures is to 
correct in-tunnei measurements at the same locations, if these are available, 
using equation (6.4). For consistency with the previous calculation 'mea- 
sured' Cp’s at the tunnel wall and ceiling 'were recreated using the source- 
source-sink model and its tunnel image system. "Empty tunnel" Au/U values 
were computed at the orifice locations and used to correct the 'measured' 

Cp’s via equation 6.4. The results of these calculations are given by the 
Pius symbols in Figures 6.10 and 6.11. The suctions are systematically lower 
than were given by the 'direct' calculation. This difference arises because, 
m applying the 'direct' method, the source and sink strengchs were not 

adjusted to an equivalent, free air value. Such an adjustment is implicit in 
equa t ion (6.4). 

There remains a question concerning the differences between suctions 
deduced from small tunnel data and those measured in the large tunnels. In 
the discussion on body shape, in Section 5, it was noted that the equivalent 
body in the smal' tunnel was slightly more slender than the free-air body, 
due to tunnel induced bubble distortion at subcritical Reynolds numbers. This 
may be the reason for the small tunnel/large tunnel differences observed in 
Figure 6.10. With the exception of one point, the supercritical data 
(Figure 6.11) agree quite well. 

It may be concluded that pressures measured at off-body points may be 
corrected via equation 6.4. If measured pressures are unavailable, the 
'direct' method may be used to generate new data from the source-source-sink 
model derived from wall pressures. The latter method is inapplicable in the 

near field because the source-source-sink body is too crude a representation 
of the actual flow. 
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7.0 CORRECTIONS TO DRAG 


7. 1 Introduce ion 

If detailed n ^osurements are available of pressures or of forces on com- 
ponents, correction on a point-by-point or individual component basis is to be 
preferred. Failing this, local information can be considered when estimating 
effective model position. However, i-^ many cases only balance data will be 
available and additional reliance must be placed upon v^/a 1 I pressure signature 
information. The possibilities in this regard are extensive but have not been 
explored fully. However, it will be seen below that gradient effects due to 
wake blockage can be estimated accurately solely from wall pressure 
measurements. 

Figure 7-1 shows the sequence in which drag corrections are made. The 
major corrections are for "horizontal buoyancy" effects which arise from wake 
blockage and to dynamic pressure for changes caused by both v,^ake and solid 
blockage. In keeping with the above sequence, subsections 7.2 and 7.3 will 
concern horizontal buoyancy and dynamic pressure corrections respect ively. 
Applications to flat plate and sphere data will be shown in subsection l.k. 
Treatment of lifting models is deferred to Section 8. 

During the present research, it was noticed that one form of the present 
correction method is equivalent to Haskell's wake blockage method (Reference 
10). This relationship will be discussed in subsection 7.5. 


7.2 "Horizontal Buoyancy" Correction 
Convert i iona I Me thods 

Empty tunnel or mode I - i nduced pressure gradient is frequently assumed to 
be constant over the model length. Buoyancy drag is then determined as the 
product of normalized pressure gradient and model volume. If the pressure 
gradient is not constant the same result may be applied to an element thickness 
6x and area A(x)..'ix to give an incremental buoyancy drag 


30 


the whole model 


6C. = - 

D S dx 


dx 


This may be written 


AC, 


A(-^ 


. dC 
X) p 

dx 


Mode 1 
Length 


dx 


AC. 


AC^(x)d 


A(x)l 


(7.1) 


-here the inlegrai includes both forward and aft facind surfaces, ac (r) the 
blockage- induced pressure on the ™odel surface, cay be determined by L of 
the methods described in Section 6 . 

'.'he "ctUQ" Method 

If model surface pressure me„ur,„„„ts are unayailable or are inadequate 
or apply, „g ( 7 .,) accurately, an alternatiye approach may be taken. As was 
seen ,n Section 5, each source-source-s ink combination determined from wall 
pressures may be used to generate an effectiye body with a defined area dis- 
tr.bu.ion and with defined pressures on its surface. principle, the body 

shape, ,ts pressures and its tunnel induced incre.»nts could be calculated 
and equation ( 7 .,) could be applied. However the Sam. result is obtained much 
more eas.l, by considering the drag of the sources and sinks themselves. 

Using the notation shown in Figure 7 . 2 , this gives 


AD 


MODEL 


= p(-Au Q - Au 0 + Au‘ Q ) 

s s w s 


(7.2) 


The equivalence between equation (7.2) and the integral of the outer surface 
pressures may be proved by constructing a control volume with an indented con- 
tour extending forward from the rear of the body and enclosing each point 
source with a sphere which isolates the singularity. Consideration of 
horizontal momentum then yields the above result. 
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If we consider the entire, closed body in Figure 7.2 and recognize that 
its (potential flow) drag is zero, we nnay write 


AD,,__ = o(-Au'^ Q - Au"^ 0 + Au Q + Au 0 ) = 

NET s s w >/ s s vj 


( 7 . 3 ) 


Equation ( 7 - 3 ) be used to simplify (7-2), giving 


AO = - pAu 


w 


We also note that Au^ is half the astmptotic velocity Q^/8H, thus 


0.^ 

AD = - ip— 


or 


AC 


0 


BH 

$ 



U BH 


(7.^) 


For non- rectangu I a r tunnels, BH is replaced by tunnel area in the above 
equat Ions . 

Applications to a Sphere 

(n what follows, a test example is set up and is used to compare various 
methods for calculating horizontal buoyancy drag. The example chosen concerns 
a hypothetical sphere over which attached flow is enforced in the presence of 
an externally imposed pressure gradient. This is induced by a set of wind 
tunnel images chosen from a typical experiment and has the general form 
illustrated in Figure 7.2. Because the flow is attached, a source-panel 
method may be used to determine the external pressure distribution, which Is 
then integrated to yield a baseline value of buoyancy drag. Four additional 
methods are eva’uated against this baseline. 
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ii) Baseline (Panel Method) 



a) Panel sphere with source sheets 

b) Apply tunnel - induced supervelocities 

c) Enforce tangency at sphere surface 

d) Determine pressures 

e) Integrate drag componenr: of external 
pressure 

Result: AC^ = 0.0216 


(ii) Conventional Buoyanoij Calculation (Figure 7,3) 


AC 


VOLUME ^ 


D REF. AREA dx 


where AC - 2Au/U 
P 


Result: AC^^ = 0.0^h0 


Hi) Use of Model Area Dis.'^Hution i Figure 7,4) 


"^0 = 


: AC (x) 

P I > 


where AC = 2Au/U 
P 


Resi 1 t : AC^ = 0.0136 at 45 pos i tion 


(iv> Area Distribution plus C 's Using Sqn. (6.4) (Hgure 7.5) 


‘S ■ - 


AC (x) d 
P 


[aU)1 

r s 


where, from (6.4), AC = C - C = (C - l) 

^ 


1 - 

* fell 

2* 


1 UJ, 

. 


1 -f 


uj 


Result: AC, = 0.0202 at 45 position 
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(v) The ''qUQ" Method 



The above examples show that conventional blockage gradient corrections 
must be applied carefully, to obtain accurate results. In particular, equat ion 
(6.4) should be used to determine surface pressure increments, since lineari- 
zation degrades accuracy seriously. Where model surface pressure measure- 
ments are unavailable, the above comparisons show that the application of 
'•pUOJ'-type methods give results which agree with the exact Douglas Neumann 
approach. Equation (7.4) may therefore be employed with confidence for 
evaluating buoyancy drag. 


7-3 Dynamic Pressure Correction 


For non-Iiftmg models, dynamic pressure corrections are applicable 
solely to viscosity-dependent drag, including base drag.' Viscous forces are 
distributed along the model length non-un i forml y and if solely balance data is 
available it is necessary to find an effective model location at which 
blockage velocity is evaluated. This may be difficult because the choice of 
location involves conflicting requirements. 
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For s,„ple configurations It may be possible to calculate effective 
»» el position with some confidence. Thus for a streamwise flat plate with 
laminar flow, shin friction varies as y-O.S ,„d it can be shown that the local 
Skin friction coefficient at 25« chord equals the value for the whole plate 

For a wholly turbulent plate, with x-0.2 dependence, the corresponding 
location 15 at 32.8^ chord. 

It is possible. In principle at least, to extend the above approach to 
more complex i„del shapes. However, skin friction and wetted periphery esti- 
mates would be needed, as Functions of x, and such an approach would be im- 
pra^ica, in „st test situations. In many cases a specific made! location is 
e by lift center, and this usually will be of overriding importance 

Thus, the quarter chord is the obvious effective iocation for an unstaPed 
simple wing. As stall occurs, it would be logical to move the effective 
ocation rearward in consonance with the lift, A constant, nominal model 
position has been retained in the cases treated so far. 

The following equation is recommended For correcting drag coefficient 
It includes both buoyancy and dynamic pressure terms 


AC. 




(7.6) 


where AC,, the buoyancy increment is given by (7.1,) and c is the value of 

Au/U^ at the model position. 


l.h Correlations for Fiat Plates and Spheres 
:^omal Flat Plates 

Figure 7.6 shows measured and corrected drag coefficients for 3 - k- 5 - 
end 6 -foot square plates tested in the XW x liV „i„d tunnel. These 'exper i- 
mental results were previously reported in Reference I. An early correction 
procedure, known as the I'emax methodi, was reported in this reference. This is 
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based upon maximum observed interference velocities. Additional corrected 
results using the "pUa- method [Equation (7.MJ are included in Figure 7.6. 

It IS evident that the newer method yields corrections which are slightly 
smaller than before. Nonetheless the corrected drag coefficients agree 
quite well with accepted flat plate drag values, though they are about .03 
higher than the value quoted by Hoerner. 

Figure 7.7 presents new results for f loor-„„.,„ted rectanguUr plates. 

Up to 0.1 S/C, the general trends are similar to before. Drag coefficients 
ere aga.n slightly higher than the Hoerner vaiue. For «4, the largest plate 
the uncorrected drag coefficient is almost double the corrected vaiue and 
the Fact that the variation with S/C is no longer linear Is probably signif- 
icant. The results show that the method can be applied successfully 

to models with cross sections of at least lOS of tunnel area and which have 
hignly separated flow. 

Spheres 

Maxv«rthy (Reference I) reviews sphere tests by a number of workers and 
shows wide variations in the dreg versus Reynolds number characteristics. In 
certain cases, it is claimed that tunnel constraint is the cause. The present 
large tunnel tests, for which blockage effects were negilgible. indicate that 
the observed venations arise From bistable flow characteristics which are 

inherent to sphere aerodynamics. Blockage-corrected small tunnel tests 
confirmed this. 


In the present tests the spheres were smooth and no attempt was made to 
tr.p the boundary layer. The turbulence level in both the small ( 30 - x A 3 -) 
and the large (l6i* x 23i') tunnels was in the O .3 to 0.5^ range, increasing 
w.th Reynolds number. The turbulence scale was smaller in the smaller tunnel 
At the t.me of testing, the spheres were consioered >'d i f f i cu 1 f models because 
repeatability was difficult to achieve. Subsequent detailed analysis showed 
bistable flow conditions over most of the test range of Reynolds numbers. In 
the small tunnel this caused speed control problems which made it very 
difficult to define transition for the large sphere. Further surge problems 
were noted with supercritical flow fully established. 
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Figures 7.8 and 7.9 show drag variation, with Reynolds number for the P 
and is-ineh spheres. «ith the esceptlon of the broken line in figure 7 8 
which shows the uncorreoted characteristic, ali small tunnel data have been 
corrected using eguation (7.6). early analyse, data appeared unrepeatable 
unre labie and generally confusing. Only after correcting the small tunnel 
ata properly, only after considering both tunnels and both spheres simul- 
taneously and only after considering point seguence as well as position did 
a clear picture emerge. The major problem in data analysis was that 
ysteresis appeared over the entire test range of Reynolds numbers. In ce-- 
tain instances -notably the ,2-inch sphere in the large tunnel - it did not 
appear because it evidently was not provoked. the small tunnel, transi- 
ona, hysteresis couid no, be detected because of speed control problems, 
he table below summarizes the hysteresis observations: 


TUNNEL 

SPHERE 

SUBCRIT. R 

e 

TRANS 1 T 1 ON 

SUPERCRIT. R 

e 

30" X ii3" 

12'' 

H? 

SPEED 

H 

15" 

H 

CONTROL 

PROBLEMS 

SPARSE DATA 

I6i' X 23 il-' 

12" 

HYSTERESIS 
NOT INDICATED 

H 

HYSTERESIS 
NOT INDICATED 



15" 

H 

H 

H 


H: HYSTERESIS PRESENT 


Against the above background, the difference, between tunnels due to 
constraint, due to tonne i - induced changes in the viscous flow or due to 
turbulence scale are minor. The 12-inch sphere data show a somewhat differ- 
ing characteristic in the two tunnels, sober i tical I y (Figure 7.8) but the 
differences are within a credible hysteresis band (Figure 7.9). Transition 
I" e e small tunnel came halfway across the hysteresis loop observed In the 
large tunnel Drag coefficient at the lower bound of the supercritical 

bistable region in the small tunnel was slightly above the observations made 
in the large tunnel. 
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The 15-inch sphere had a slightly more 'peaky' subcritical drag char- 
acteristic than did the 12-inch sphere but large tunnel and corrected small 
tunnel data defined the same hysteresis boundaries in this region. Transi- 
tional and supercritical data were sparse in the small tunnel. There was 
reasonable agreement regarding transition Reynolds number but large tunnel 
supercritical drag levels appear high. The small tunnel supercritical 15 " 
sphere drag level, after correction, agrees with the large tunnel 12" 
sphere data. 

Because of the difficulties just described, it is tempting to suggest, 
in retrospect, that the sphere was not a good candidate for checking the 
blockage correction procedure. However this is not so, ' ecause the diffi- 
culties encountered here would also be found in many wing-stall situations, 
for example. The preceding analysis shows that the new correction method 
worked well and was, in fact, the least of the problems encountered. Cor- 
rected data from even the large, 15 -inch sphere in the 30" x 43" tunnel is 
considered at least as reliable as the large tunnel data at subcritical, 
transitional and at supercritical Reynolds numbers. 


7-5 The "Crnax" ^nd Maskell Methods 

Reference 3 describes two early methods of implementing equation 7.6 
which have now been superseded. One of these, the method, gives very 

similar results to the present, recommended method but is less soundly based. 
Using the "cmax" procedure, the peak interference e„,ax replaces e^. The 
buoyancy term is determined from model volume and pressure gradient, rather 
than via the “pUQ." method. 


A velocity vector Vq(x) is accurately measured in the wind tunnel and 
uncorrected pressure coefficient Is determined via 


C 


V (x) 
o 


U 

CO 

u 


2 


(7.7) 


“eW - ^ 0 * .(.» 

Equation (7.8) reorecani-c 

Equation ( 6 . 4 ) The "e - '""T P''®ssure correction given by 

W procedure is not x-depe„de„t end (7.8) reduces to 


U (x) = u (l + e ] 
® ^ maxJ 


( 7 . 9 ) 


The corrected pressure coefficient r t, • . 
effective „einstree» speed, is ,il 


fvjx)' 

(irra-i 


Eliminating v^(x) between (7. 


7 ) and (7. 10) 


g I ves 


( 7 . 10) 


i<J„{x) 


in Maskeil's notation. 


!c _ ' ' 'p„ 


( 7 . 11 ) 


T-his equation, which derives from the "e a • • 

Haskell's equation (17). The -c '• ru a 

Haskell's method but replace! but^e " 

assumptions, by results derived ^-m Iq^^rinrwT^ 
measu^erT.Jn^s, p^ehsure signature 


epy change in bubble shape (wh ic‘h“I''hevnTw''L' 

which base oressnr® ^ small; . For cases in 

uabe pressure is not measured Maclddaii , 

bubble area change is (B/C) f i ! 

•n9c .s (B/C) t,™es that for a full, vieiding wake, where B is 
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bubble cross sectional area and C is tunnel area. This changes the relation- 
ship For determining Haskell’s k values but is automatically accomrr,odated by 
the wall signature procedure. The present method implicitly includes bubble 
distortion effects while avoiding Haskell's possibly contentious 
assumptions in this regard. 


0 CORRECTIONS FOR LIFTING AND POWERED MODELS 


8.1 Flat Plates at Angle of Attack 

Linear- - a Range 

Figure 8.1 shows corrected lift curves for four flat plate wings tested 

o just beyond stall. The semi-span models were mounted from the tunnel 

r oor and the whole-wing aspect ratio was 3 - Plate area ranged from 1.6| to 

^o.//o of test section area and the "e " nK=.f-hr^^ , j 

une method was used to correct 'q' 

standard Kxt book «,hoda war, used for angl, of attack corrections. These 

were senerally less than one degree. 

Despite the very Urge «,dei site for wing Wk (circles) and despite the 

presence of an uncontrolled boundary layer on the tunnel floor, the lift cor- 

relat.on between the four oodeis is excellent over the linear C, - „ range 

here was so«e tendency for smaller oodels to stall earlier, which was prob- 

a y ue to the relatively greater adverse influence of the floor boundary 

aver. The li.t curve slope is about 10* greater than that given by simple 

wing theory; this is most likely a manifestation of lift from edge vortices 
at the wing tip. 

the f, "'7° ID this case the effect of 

h oor boundary layer was mare serious. The percentage effect is smallest 

Urge wing, yi, and C„,/Ctt for this wing agrees well with the free air 
-iue guoted by hoerner (broken line). The ,/gii term in Hoern.rU eguation 
represents the loss of leading edge thrust for the thin flat plate The 
smaller wings exhibit a ..cut oiif effect in the root region due to the tunnel 
oor oundary layer. Hoerner states that, with 10* of span removed centrally 

the M " Even for 

s^ est wing tested ,5,5.. semi-span), the effect is significantly Uss 


The above results show that the combined blockage and conventional 
angle of attack corrections produce good results for large wings. Including 


loor boundary layer, 


one for which Cl S/C was 0.13^. In the presence of a f 
the largest semi-span model gave the best results. 


Pos't—CDtaL'L Rcch'^g 


.gures 8.3 and 8.4 extend the previous lift and drag data all the way 
to 90-degrees angle of attack, i.e. the normal flat plate case. Uncorrected 
are also shown. Drag corrections are particularly large at the higher 
angles of attack. The correlations between models remain good except in the 
immediate post-stall region. Here, it appears that the floor boundary layer 

again had an adverse effect on the smaller models, this time in delaying the 
approach to the second lift peak. 


F.gura 8.5 shows the above deta io polar for« at a larper scale. This 
raes.phas.tes the large wagoitude of the blockage corrections at h.gh angles of 
atwck. It ,s also apparent that the method overcorrects slightly in 

th.s region for the largest wing. On reverting to the "pUO" method (tagged 
cricles), improved correlation is obtained up to about 60-degrees angle of 
attack but thereafter results from the "bUQ" and methods converge. 

Examination of derived blockage parameters reveals that model-c, e„. and 
wake source strength all Increase smoothly with angle of attack in this range. 
However, the peak blockage position moves back up to 60-degrees angle of 
attack and then stays constant. This evidently is what causes corrections 
derived by the "pUQ" and methods to converge. 


The above comparison between the two correction methods 


portant distinction related to lift. The e 


raises an im- 


max 


method, by imp) I cat ion, 


includes a horizontal buoyancy term in application to drag which appears as 
the diiference between and model e. This same difference is also applied 

to hft. This is incorrect because horizontal buoyancy cannot induce lift. 
Th.s ,s a further reason for favoring the "pUQ.. method over the method. 

Figure 8.6 shows pitching moment data over a 0 to 90° angle of attack 
range. The correlation between the four wings is similar to that for lift and 
drag. At 90-degrees the quarter chord pitching moment equals 0.25 times the 
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drag coefficient. Observed differences here reflect both drag coefficient 

variations, between models, and possible small set-up errors for the smaller 
mode 1 s * 


8.2 Correlations for the Knee-Blown Flap (KBF) Model 


in its basic form, with slats added but no tips, the reference area for 
the KBF model is 6.2^ of the 30- x tunnel cross section, or 37% of the larg 
est flat plate wing. However, the Cj_-range is an order of magnitude greater, 
so C|_-S/C and hence angle of attack corrections are about three times greater 
than for the fiat plate wing. Figure li.7 of Reference 2 shows that angle of 
attack and blockage corrections are of comparable magnitude for the KBF model 
in the 30" x 43" tunnel. The angle of attack corrections used are those of 
Williams and Butler (Ref. 12) which are applied as described in Reference 2. 


Baseline data were obtained in the NASA/AAMRDL wind tunnel (see Refer- 
ence 2). Since angle of attack corrections are significant, comparisons 
between corrected small tunnel and large tunnel data are as much a test of the 
angle of attack corrections as they are of the blockage method. There is a 
need to carry out both corrections on a consistent, unified basis. This is 
discussed further in Section 9. 


Wall pressure signatures reflect not only model flow conditions but 
also the effects of changes in tunnel conditions, particularly flow breakdown. 
This is reviewed in some detail in Section 6 of Reference 2 which discusses 
signatures for moving ground, blown BLC ground and very high cases. 


Figure 8.7 shows the boundary for jet impingement on the tunnel floor 
as a function of Cjj and a. A firm correlation has been established between 
the occurrence of impingement and the subsecunnt development of a second peak 
in the wail pressure signature. The broken line in Figure 8.7 is the boundary 
beyond which the seco.nd peak dominates the pressure signature. Once this has 
happened, the correction process is seriously compromised and model flow con- 
ditions are probably unrepresentative of free flight. Figure 8.8 shows the 
lift characteristics of the KBF model across the entire test range, to C ' s 
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curves are 


of 10 and a-s of 30. Individual points are not shown because the 
the product of cross plots needed to obtain constant values from corrected 
data. Up to Cp - ,, there is generally good agreement between the two tunnels 
except that unpowered lift is slightly higher in the small tunnel. 

Corrections are very small in this case. 

The C, = 2 curve is interesting because it intersects both the impinge- 
ment and the strong second peak loci in Figure 8 . 7 . The large tunnel and the 
corrected small tunnel results start to diverge half way between the two con- 
ditions. At Cy's of k.O and greater, differences between the small and large 
tunnel data increase as impingement intensifies. 

Figures 8.9 and 8.10 show drag and pitching moment correlations corre- 
sponding to a Cl range up to 10. Drag correlations are nuite successful, 
including the unpowered case, but overcorrection again occurs at C =4, the 
impingement case. The trends in the pitching moment correlations (Figure 
8. lO) are generally similar to those for lift. 


Reference 2 includes correlations similar to those described here but 
using what is, in effect, an early version of the present correction method. 

A "true Q" systenr was e.ployed which relied upon e single pressure meesured 
pert way down the rest section. Exaniination of wall pressures measured sub- 
«,uentiy has shown that the orifice used in the "true system lies between 
she rirst and second peaks of this signature under most circumstances, which 
was a fortunate choice, in effect, the previous system isolated the wake 
blockage from the total signature and ignored the solid blockage. For this 
reason the new procedure gives better results up to and including c„ - 2. 
Beyond this, the undercorrection tendency of the old, "true q", procedure 
lead to somewhat better corre 1 a t i ons , 


Some limited studies have been made 
with the second peak in impingement cases, 
due to a second effective body situated at 
location, this body retards the Flow In the 
signature was calculated for a Cjj = 8 data 


of blockage interference associated 
This peak may be interpreted as 
floor level. Because of it‘s aft 
region of the model. A residual 
set by subtracting the fitted, 


peak SKjnacure fro™ the ™easered one. A toosh analysis of this n.sl- 
oe 9 ave brockajie increments which improved the lift correlation Ho 

9- analysis provides a likely explanation for small tunnel effects and 
m.9 . e regarded as a basis for further correction, any results obtained 

speculative, at best, because of impingement- I nduced flow 

d i s tort ion . 


9.0 DISCUSSION 


With methods for determining axial velocity increments due to blockage 
now weil established, emphasis shifts to determining their effects on a self- 
consistent basis. need has also been noted for tunnel surface pressure 
signature analysis method which determines corrections to angle of attack. 

To achieve this, a simple horseshoe vortex model might be used, which com- 
plements the line source-sink model employed for blockage estimation. The 
method could be implemented via look-up charts with inputs from tunnel 
ceiling and floor pressure signatures. Such a procedure would define the 
angle of attack environment of the model, including variations in at least 
the axial direction. The problem of determining the effects of the 
distributed angle of attack increments parallels that concerning distributed 
blockage velocity increments. 

Significant advantages accrue when blockage and angle of attack con- 
straint corrections are determined on the above unified basis. For example, 
blockage increments found via the present method could be applied to the 
horseshoe vortex model determined during the angle of attack analysis. Tli i s 
would permit lift and possibly pitching moment increments to be determined 
using information derived solely from tunnel surface pressures. The present 
need to estimate a proper effective model position would thereby be avoided. 

It IS implied in the above discussion that the vortex lifting model 
would comprise one or maybe two horseshoe systems of unknown strength 
situated at unknown positions, i.e. a system having non-linear geometric 
unknowns and linear strength variables as for the present blockage analysis. 
This approach would probably work well for small models or straight wings but 
for more complex configurations more elements may be needed to attain the 
necessary resolution when calculating induced effects. A preferable approach 
might employ two linear systems, one each for the source and vortex flow 
models, each of which has multiple elements at known positions but of unknown 
strength. "One shot" analyses would then be possible using predetermined 
influence matrices and iterative techniques would be avoided. 
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The above discussion is of course, quite speculative and any attempt to 
implement the suggested scheme should probably be made on a careful, stage- 
by-stage basis. However, much of the data described in the present report 
could be used for checkout purposes. 


10.0 CONCLUDING REMARKS 


The present report concerns a wind tunnel blockage correction method, 
first described in Reference 1, which relies upon wall pressure measurements. 
This ensures that changes in blockage, due to stall or power effects for 
example, are sensed immediately and are incorporated au toma t i ca 1 1 y into the 
correction process. Correlations from tests on various models are given in 
this report including normal flat plates, spheres, an idealized automobile, 
flat plate wings and a blown-flap model. Tests upon scaled models in a 
single tunnel or upon one model in small and large tunnels are included in 
a 1 1 cases . 

The terms of reference for the above studies relate predomi nant I y to 
the prediction of b 1 ockage- i nduced velocities within the test section volume 
Secondary objectives concern the conversion of such velocity increments to 
model surface pressure or force corrections. 

Determination of Interference Velocities (Section 4) 

The development of wall pressure analysis techniques and blockage 
velocity prediction methods is complete. Velocity increments due to combined 
wake and solid blockage have been computed throughout the test section volume 
and it has been found that knowledge of the centerline axial velocity incre- 
ment is sufficient ^or most purposes. Certain ’“problem" signatures —notably 
those with double peaks due to tunnel flow breakdown —are accepted and 
analyzed by the current program. However the validity of the test itself 
should be questioned when such signatures occur. 

Corrections to Measured Fressia^es (Section 6) 

A correction technique has been devised fcr model surface pressures 
which relies upon a "local mainstream" concept (see equation 6.^). The 
method has been checked out successfully against definitive theoretical pre- 
dictions, in tunne 1 -versus- free-a i r sphere surface pressure comparisons, and 
in automobile tests* It is suggested that a similar concept may also apply 
to. skin frict ion. 


^8 


^cr>veotions to Drag (Section 7) 


Dra, correction t.ch„i<,„os have been r.vie-ed in so»e detai , An earlv 

technique, designated in Reference I as the "e " merhor) • i ^ 

to Haskell's hirAv-iv lethod, is closely related 

blockage correction method for bluff bodies and stalled „ings 

"ore recentlv. the -o«<r .ethod has been developed which involves separaL' 

huovancv and dynaoic pressure corrections (see eguation 7.6). euovancy 

tenn.ned at the ™odel center. The two .ethods produce co»parabie results 

» work wen for non-liffin, .^d„s. However, it is shown that the p 

-S6od ,s better founded and therefore should be oore readilv app, i a . to 
complete, I ifting models. PP-'caoieto 


Corrections for Lifting Models (Section 3.1) 


The above methods extend to lifrinr, v-co 
nrr provided that 

"tt IS concentrated around a sinqie x-lorarinn tk 

method bacons evident in this ann, , " 

«av be d ■ ■ , application s.nce dynamic pressure corrections 

Half-Chord under stalled conditions. Tests on a family of * 3 flat plate 
■ngs up to and including 90 -degrees angle of attack showed excel, .„t corre- 

fPTWingitunnel area ratios of up to ,0?. There were up to S* drag 

discrepances above 60-degrees for a wing with S/C of .167 but the correctLs 
were extremely large - the uncorrected drag coefficient being twice 
erected value with the wing norma, to the flow. An effectL noder 

alg7„°f 7a“k -onal 

angle of attack corrections were employed. 


•^orrecttons for Powered Models (Section 3.2) 

approx "7T T'" “"'P'Pd successfully by the present method up to 

approx, mately Cg . ,0. c„ - A, despite the fact that mild tunnel flow break- 
oow„i„as present for all c„ - A data. The occurrence of flow breakdown Is 
y etecte from the wall pressure signatures which developed a second 
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peak under these conditions. This caused increasing ove rcor rect i on , by the 
present method, as increased beyond 4- 

Modei Size Const dem.t ions 

In the past, model size in relation to the tunnel cross section has 
sometimes been limited by our ability to perform blockage corrections with 
confidence. Development of the wall signature method has removed much of 
this difficulty and much larger models may be considered. An extreme example 
is described in Section 7 concerning a 15”inch diameter sphere in a 30" x 43" 
test section. Both surface pressures and drag values were corrected 
successfu I ly . 

Axial gradient of blockage velocity may be a limiting factor for more 
complex models, particularly if lift is present. If lift is not concentrated 
around a single x^location as on a large^chord wing, a highly swept wing or 
a whole aircraft model with a tail — accurate moment corrections may not be 
possible unless pressure d i s t r i bu t i ons or component forces are available. 

Thus b 1 ockac e gradient, rather than bl ockage magn i t ude , may 1 i m i t nK)de 1 size 
in certain cases. Because of properties inherent to the present method, 
gradient presents less difficulty regarding axial force corrections than for 
transverse forces or moments. A unified correction method is proposed, in 
Section 9, wnich extends the present v>/al] pressure signature analysis method 
by adding .c-.ex elements to the source-sink theoretical model, thereby 
placing ar.j’i of attack corrections on a similar footing to the present 
blockage cm 'ection. 
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Figure 3.1 Summary of blockage tests. 













Figure 3-2 Details of normal plates tested in the Lockheed 
t6i' X 23 i ' tunnel . 
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Figure 3-3 
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Figure '».! Ceiling pressure signatures for a floor-mounted 
semispan wing (Fig. 3 . 3 ) at 0°<ai90® (S/C = 16.?^) 
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(a) Measured Ceiling Supervelocity 
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Buildup of tunnel centerline interference velocitv 
from individual source and sink contributions. 
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EFFECTIVE BODIES FOR A SUBCRITICAL SPHERE 
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Figure /.8 Ldrcjc- Lunne 1 and corrected sfiia 1 i - t unne I drag coefficientb 
for a 12-uich sphere. 
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Figure 8.10 Large tunnel and corrected small tunnel pitching 
moments for knee-blown flap model. 
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PROGRAM DESCRIPTION 


APPENDIX I 


PROGRAM DESCRIPTION 


The blockage program was initially developed for the Lockheed l6^-' x 
23i‘ low speed wind tunnel using, initially, its CDC 1700 computer and later 
its Tl 990 system. Complex pressure signatures, arising from powered model 
testing, in the Lockheed 30“ x 43“ tunnel led to further refining of the 
program, using that facility's Lockheed MAC-16 computer. The program has 
been further tidied up for implementation on the NASA CDC 7600 system. 


This appendix gives an overview of the program together with the 
input and output information needed for its use. No reference to 

the main text of this report should be needed in order to run the program. 

Listings of the program are given in Appendices II and III, and the required 
working charts may be found in Appendices IV and V. The Information avail- 
able on the 7600 consists of the program itself and three data packages 
(working charts for the appropriate tunnel configurations). This information 
is stored under the identification CORSIGLIA and the following names: 
NAS 2 - 9883 A Program and library. 

NAS 2 - 9883 B Charts for .699 aspect ratio rectangular tunnel, 

NA$2'-9883C Charts for 40' x 80' center tunnel inount. 

NAS 2 - 9883 O Charts for 40' x 80 ‘ semispan floor mount. 


General Structure . - The general structure of the computer program is 
shown in Figure A1*, A listing is given in Appendix M. The program can be 
conveniently considered in three parts: input, computation and output. A 

fourth part, not shown, is the library. This consists of miscellaneous 
routines, some of which are standard systems library features on the 7600 . 
They were required for the Lockheed m i n i -compute r , and have been retained in 
anticipation of the program eventually being transferred back to a mini- 
computer. A listing of these routines is given in Appendix IN. 
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Basically, the program accepts the tunnel wall pressure signature and 
subtracts from it the empty tunnel signature, leaving the signature due to 
the presence of the model. This signature, defined by a pressure coefficient 
distribution, is converted to a velocity increment distribution. 

The velocity distribution is split into a symmetric and an anti- 
symmetric component. The source-sink pair that would produce the symmetric 
component and the source that would produce the antisymmetric component, are 
derived from the included working charts. It is then possible to calculate 
the velocity distribution, along the tunnel centerline, that results from 
tunnel effects stemming from this source and source-sink pair combination. 

The resulting blockage parameters are displayed. 

The velocity distribution along the tunnel wall, also calculated from 
the combination, is displayed graphically, on the line printer, together with 
the test data to show the degree to which the program is matching the test 
data. 


Two add i t i on a 1, optional, printouts are available. The first is output 
from within the computational portion of the program to facilitate debugging. 
This may also provide a better understanding of the program's attempts to han- 
dle a non-standard velocity signature. The second option is a breakdown of the 
the components that make up the wall and center tunnel velocity distributions. 

J_np^t_. - The required input to the program can be considered in three 
parts. These are, initially the working ciiarts, secondly the main input and 
finally the run input. 

The input format for the working charts can be seen in the program list- 
ing (Appendix II) in the subroutine EPSTAB. The data blocks provided (named 
in the beginning of the appendix) are complete and in the necessary format for 
immediate insertion into the program. 

Subsequent input is shown in the program listing in the subroutine 
EPSINP. 


Mai^'Z Input 


The main input portion of this consists of the following cards: 

(1) IPR ITEST 215 

I PR - Output device number (line printer) 

ITEST - Test number. 


(2) ATUN BTUN 2F8 

ATUN - Tunnel cross-sectional area (ft-) 

BTUN - Tunnel reference width (ft) 

For the Ames kO'- x 30 ' tunnel, BTUN = 80 and ATUN = 2857 for 
center tunnel tests and ATUN = 5713 for floor-mounted half 
mode I s . 


(3) TITLE 

TITLE - Test description 


20A4 


ih) 


BMOD 

BFAC 

SHOD - Model 

span 

8 FAC - Ratio 

of SOL 

XMOD - Axial 

pos i t i 


XMOD 3F8 

(ft) 

irce length to model 
on of model (ft) 


span (assume 1 . O) 


(5) NWST 15 

NWST - Number of pressure readings (maximum = 21) 


(6) XWST(1) XWST(2) XWST(NWST) iOF8 

XWST - Array of pressure tap axial positions (ft) 


1 I 1 


(7) 


CPEM( 1) 

CPEM(2) 

I 

! 

CPEM(HWST) 


CPEM 


10F8 


Array of empty tunnel pressure coefficients. 


( 8 ) 


XMC 

XMA 

1 

VOLM 1 SF 

AR 

DELTA 


6f8 


XMF - Forward extremity of model location (ft) 

XMA - Aft extremity of model location (ft ) 

VOLM - Mode! volume (cu. ft) 

SF - Wing area (sq. ft) 

AR - aspect ratio | 

r for incidence correct ons 

DELTA- interference factor 

The interference factor ( 5 ) is that obtained from Pope (ref. 7) 


(9) 


LU 

lUSES 

lUSEE 

ILIST 


IDEBUG 


I TAB 


615 


LU - Plotter device number (LU=0) 
lUSES - Forward asymptote identification 
lUSEE ■■ Aft asymptote identification 
ILIST - Additional output option (1L1ST?^0) 

IDEBUG - Debugging output option (iDEBUGj'O) 

ITAB - Working charts output option (1TAB?^0) 

The testing of floor-mounted models requires that the image sysLcr.i (the 
floor being the axis of symmetry) be recogni.red. Tunnel area, refer- 
ence width, model geometry and interference factor should include this 
I mage . 


LKrut 


(I) 


ALPHA 


F8, F6 


ALPHA - Angle of attack (degrees) 


I 12 






(2) 

CPWST ( 1 ) 

CPWST(2) 


I CPWST (NWST) 

CPWST - Array of pressure 

' j 

coe f f i c 1 en ts 


I0F8 


(3) 



(JU - Uncorrected tunnel 
PU - Uncorrected tunnel 


2F8 

dynontic pressure (psf) 
static pressure (psfa) 


i^) 

CODE 

CMUU 

CLU 

CDU 

j~ CMU 

CODE - (Run number) + (Point number) /lOO 



CMUU - Unco rrec ted power coefficient 


CLU 

Uncorrected 

1 i n 

coe f f 1 c i en t 

CDU 

- Uncorrected 

drag 

coef f i c i en t 

CMU 

Uncorrected 

pitching moment coefficient. 


5f8 


(5) 


lOPT 


15 


lOPT - Next input option (-1 Terminate. 0 New run input. 


A sample Input is shown in Table 1 


NevA/ norn input) 


Computation . - An outline of the computation procedure is shown in 
Frgure A2 . Practical and theoretical considerations of the measured velocity 
distribution show that a similar distribution can be obtained from a distri- 
bution of line sinks and sources along the centerline of the tunnel. The 
purpose of this section is to determine the source and source-sink pair that 
would generate a signature matching that measured. The effect of the source 
and of the source-sink pair is to generate two velocity distributions which 
are respectively antisymmetric and symmetric about a plane in the 
approximate vicinity of the model. 
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The amplitude of the antisymmetric portion (UOUMIN) is determined from 
e orward and trailing asymptotes of the measured velocity distribution 
)• The axis of antisymmetry (XVOB) is initially assumed to be at the 
oael position tXMOS) and the distribution determined (UAOU) . Subtraction of 
e istribution from the test data leaves the symmetric component (USOU). 

A parabola is fitted to the upper portion of the symmetric data, using 
east squares fit routine, and the axis of symmetry (XPOB) determined. If 
's IS appreciably different from the assumed position (XVOB) the anti- 
symmetric distribution is recalculated, at the new position, and the above 
process reDeat:ed. 


The source (QR-/) that would generate the antisymmetric distribution is 
c.rectly oroppional to the amplitude of this distribution and is positioned 
axis o antisymmetry. The source-sink pair, strength (QFS) and 
spacing CSOB , are determined from the working (or look-up) charts and are a 
unction of the amplitude of the symmetric distribution (UOUMAX) and the width 
of the distribution at the half-amplitude height (DXOB) . 

This information is now sufficient to calculate the theoretical distri- 
bution of velocity at the wall (UOUW) and along the tunnel centerline (SIGUOU) 
Not Shown on the How chart of Figure A2 is a check on the validity of the 
uied asymptotes. If conditions are not asymptotic at the selected forward 
.nd aft d.ta poinls (dpun„i„ed b, .USES and lUSEE), then the calcdiated 
curve (UOIJW) will not pass through these selected data points. This is 
checked at the end of the first pass through the calculation and the 
asymptotes modified by the amount of miss. 

Provision is made for the program to byp.ss the source-sink estimation 
■n the event that it is unable to adeguateiy curve fit the USOU data ard 
proceed to a subroutine named PUNT. In this case the blockage, a, the model 
position, i, estimated directly using a yersion of Hensel's method (ref. 5) 
and no distributions are available. 


I U 




Additional information on the program is available in the form of 
comments on the listing in Appendix II. 


Debug Output . - The setting of I DEBUG^O will result in output similar 
to that shown in Table 2. The first line (U/U) is the test data. The curve 
extremities are to be defined by the first point (IIJSES=l) and the third from 
last point (I USEE=-3) • 

The next six lines represent one pass through the iteration loop. A5~7 
represents the three coefficients of the antisymmetric curve A5, A6 and A/. 
UA/U and US/U are the antisymmetric distribution and the residue (U/U - UA/U) , 
The next line, labelled XV/B gives in turn, the assumed axis of antisymmetry 
(XVOBK the coefficients A5, a 6 and A7 , the coefficients of the parabola A, B 
and C, the amplitude of the symmetric curve (UOUMAX) , the width at the half- 
amplitude (DXOB) and the location of the axis of symmetry (XPOB). 


At the conclusion of the first pass through the iteration loop, the 
program bypasses the iteration check and calculates the total wall distribu 
tion (UOUW) and checks the validity of the selected asymptotes. The itera- 
tion Output is repeated. Note that the forward asymptote (point 1 in the 
UA/U data) has been decreased by .0006 and the aft asymptote by .0001. The 
axis of antisymmetry and symmetry was initially assumed zero and calculated 
to be -.0327. On the second pass an assumption of -.0327 gave an answer of 
-.0351 which is within the required limits and the iteration terminated. 

The next line of data, labelled BS/B, gives the inputs to the working 
plots, BSOB and DXOB together with the results from Chart I and Chart II. 
The final line of data, labelled XM/B, gives the extremities of the model 
(inputs), and the velocities at these points, and finally the buoyancy 
factor (DCPDX). 


Normal Output . - An example of normal output is shown in the upper half 
of Table 3. The Symbo.s on the plot are 'O' fa-- .he test dsta and 'X' for 
the matching calculated distribution The supplied coefficients are shown as 
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measured, with blockage corrections and with incidence corrections per Pope 
(ref. 7). 

The following data consists of the blockages (EPS) at the model position 
[X(H0D)/B] and at the position of maximum blockage [X(MAX)/B], the source-sink 
strength (QS/UBM) and the single source strength (QW/UBH) . The model span 
(BS/B) and the width at half-amplitude of the symmetric velocity distribution 
(XV/B) provide the inputs for the working charts, giving the source-sink 
spacing (CS/B) and the source-sink strength term (UFM) . US(MAX)/U is the 
amplitude of the symmetric curve. A5 , A6 and A? are the coefficients of the 
antisymmetric curve and DCDWB is the buoyancy factor. 


The comment on the final line indicates the status during the iteration 
loop where divergent conditions or a jump to a Hensel computation (subroutine 
punt) may have been experienced. 


Additional Output. - The lower half of Table 3 shows the additional 
output that is optionally available. It consists of the pressure point 
locations (X/B) , the test pressure coef^cients (CP) and the resulting velo- 
city distribution (U/U) . The antisymmetric component (UA/U) , the residue 
(US/U) and the calculated wall distribution (UOUW) follow. The final three 
sets of data are the center-tunne I build-up of the potential blockage, 
viscous blockage and total blockage. 
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Table 2 Sample debugging output. 


U/U FULL SCAL6 ■ <*033? 



Table 3 Sample program output. 
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PROGRAM LISTING: AERODYNAMIC ROUTINES 
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PROGRAM LISTING: LIBRARY ROUTINES 
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WORKING CHARTS FOR RECTANGULAR WIND TUNNELS 


APPENDIX IV 


WORKING CHARTS FOR RECTANGULAR WIND TUNNELS 


The later stages of wall pressure signature analysis require references 
to data charts or carpets which are used to determine blockage line source 
and line sink details from measured signature character I st ics (see Figure 
2 . 3 ). For rectangular tunnels, the influence coefficients needed to define 
these carpets can be written .n closed form as given below. Less regularly 
shaped tunnels, such as the ' x 80 ' , may require a numerical approach (see 
Appendix V), 

Equations for a Soux^ce-Sink Flow Model (Rectangular Tunnel) 

For a flow model involving three horizontal line sources on the center 
plane of a rectangular tunnel, the axial velocity increment induced at a 
point X, y, z is giver, by 


Au = ^ 

L=2 


Sits 


MM h 

I I R W 


(x - xp 


I u u 2BH ABS (x - x7y 

L m=-M n="'N L=2 L 


h 

I 

L=2 




4 

I 

L=2 


2 BH 


(A'l) 


In equation (Al), is an axial velocity Function, for axial velocity 
at a receiving point x, y, z, given by 


R = 


X 

{x2 + z^} 


(Y + s^) 

{X2 + (Y-S|_)2 + 


{X2 + 


(Y - s^) 

(Y-s^)2 Z2}^_ 


whe re 


X = X - X. 

L 

Y = mB - y 
Z = nH - 2 . 


L is the index for the source system and m, n represent the image cell loca- 
tion in the cross flow plane (m = n = 0, representing the central source, is 
included in the summations when wall velocities are calculated but excluded 
v\/hen interference velocities are required.) L = 1 Is reserved for the model 
I oca t ion . 


J, in equation (Al) is obtained from 


J = tan 


- 1 






X{X2 + + Zi2}i 


- tan” 


X{X" + Yi- + Z2“} 
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- tan 


- 1 


Y2 Zi 


Y, Z- 


X{X^ + Y2" + Zi-]‘ 


+ tan 


X{X2 + Y.- + Z-)-}^ 


whe re 


(x - xp 


'1 = (y + (2M+1) |) 


Y^ 5 (y - (2M+I) 


Zj H (z + (2N+I) y) Z 2 = (z - (2N+I) y) 


The tangent function^ 7* concerns the velocity induced at the “receiving" 
point X, y, 2 by a transverse rectangular source plate si stated at x = x^ 
and covering an area (2M+l)B by (2N+l)H. 

For each source system, L in equation (Al), the line source function 
R, is summed over a finite nu.iber of tunnel images in the first term. The 
middle two terms approximate the region outside of this by considering it as 
a doub 1 y- i n f i n i te source sheet (second term) minus the “hole" (third term) 
which is already represented. The final term represents downstream sink 
sheets needed for continuity and ensures zero disturbance far upstream of the 
mode ) . 


U2 


As written, there are nine unknowns in equation (Ai); three source 
strengths, , three line source positions, X|^, and three source semi-spans, 
The problem is converted to more convenient notation and reduced to 
seven unknowns by putting 


and 


Q3 = +Q, "3 = 2 ^ 

Q4 = -Q 3 J 1^3 

^ ^2 = I ‘^W 


(A2) 


(A3) 


The subscript s refers to solid or separat ion-bubble blockage and w refers 
to wake blockage. 


Prepaixttion of Working diarts (a) Sign.iturr 


:gsz3 


(Charts I and II) 


Charts I and I I are used v^;hen matching the measured signature to the 
idealized mathematical model of the flov^/. Prior to using these charts, the 
measured signature is resolved into symmetric and antisymmetric parts. 
Charts I and II are then used to analyze the symme tr i c ,be 1 1 -shaped part. 
P’^eparation of these charts is described briefly below: a more detailed 

description is given in Reference 1. 


Chart I (Pjgure A3) determines source-sink spacing, c^/B, from the 
assumed source and sink spans b 5/8 and the width of the symmetric part of the 
signature at half-height, Ax/B. To prepare the chart the symmetric terms, 

L = 3i ^ lo Equation (Al) are evaluated at the wall location over a range in 
x/B between ±1.0 for the c^/B, b^/B combinations of interest. The central 
source and sink, for which m = n =* 0, are included in the summations. For 
each C 3 /B, b 5 /B combination, the peak width at half-height, Ax/B, is deter- 
mined, The results are then assembled suitably for chart or table look-up 
with Ax/B and b^/B as independent variables and Cg/B as the dependent variable. 
Cross-plotting is used to obtain "whole“ values of Ax/B for this purpose. 


1^3 


Chart I! (Figure A4) is used in the determi na t ion of source and sink 
strengths from the same combinations of input variables, 65/6 and 

Ax/B. This is a departure from the Reference 1 procedure which employed Cs/B 
and Ax/B. To prepare Chart II, the symmetric terms in Equation (Al) are 
evaluated midway between the source and sink, i .e. at the velocity maximum. 
The quantity (riniax^/Qs^ determined for each Cs/B, b^/B combination. For 
conventence and accuracy, c^/B is replaced by Ax/B for chart or table look-up 
purposes . 


Preparation of Working Charts (h) Timnel Centerline 
ryiterference fC?ta 2 ^t III) 

"Chart 111“ (Figure A5) gives the normalized tunnel centerline inter- 
ference velocity, due to a single line source, as a function of x/B. Inter- 
ference velocity is calculated using L = 1 and omitting m - n = 0 from the 
summation in Equation (Al). Figure A3 is drawn for a b^/B value of 0.30, 
Strictly, a carpet should be drawn to include other bg/B values. However, it 
is found In practice that dependence on bs/B is very weak and the curve shown 
is adequate throughout the bs/B range of interest. 

The total source-source-s ink interference signature is built up from 
rigure AS using suitable scaling and superposition procedures. 

F loo r^Mo i in ted Ha If -Mode I s 

For a rectangular tunnel of app rox i ma te 1 y vT aspect ratio, the 
geometry is the same for a floor-mounted half model as it is for a centrally- 
mounted model. The same charts and procedures may therefore be used provided 
that tunnel height and width are increased so as to represent the “double 
tunnel" properly. 



igure A3 "Chart I"; Determination of source-sink spacing 
(Rectangular tunnel, H/B= 0.699). 
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"Chart M": Determination of source-sink strengths 

(Rectangular tunnel, H/B= 0.699). 
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WORKING CHARTS FOR THE 40 ' x 80‘ WIND TUNNEL 

The image methods described in Appendix IV cannot be applied to the 
AO X 80' tunnel cross section because of its unusual shape. Errors for 
approximations to the shape (e.g. ellipses) were difficult to quantify and a 
detailed, surface panel technique was therefore employed. 


Preparation of Ctharts I and II 


Figure A6 shows the theoretical flow model used to generate data charts 
for the AO' x 80' tunnel. A section of tunnel about 3-5 diameters long is 
panelled with doublet sheets. Velocities due to a centrally located line 
source [Figure A7(a) ] are calculated at each panel center. Various line 
source spans are used in successive calculations. In order to assure correct 
asymptotic behavior far upstream and far downstream, velocities due to ring 
vortex tube extensions are also calculated and added to the source- induced 
velocities. The test section boundary conditions are now defined completely. 
In the next stage of the calculation, an influence matrix is computed for 
the doublet panels and doublet strengths are found by inverting this matrix 
and applying the boundary conditions. Tunnel wall supervelocities are then 
computed and the resulting curves are developed to give "Chart I" (Figure A8) 
and "Chart II" (Figure A9) in the manner described in Appendix IV. 


Preparation of "Chart III" 

"Chart III" (Figure AlO) gives the axial velocity along the model cen 
terline with the central source removed. For this, velocities due to t.ie 
ring vortex tubes are added to velocities induced by the doublet panels 
representing the test section. 


In all of the above calculations the bilateral symmetry of the test 
section is exploited in order to reduce matrix size. Although there are 50A 
panels on the test section, only a 126 x 126 matrix is inverted. 
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Ploor^-Mointed Half Models 
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